16.451 Lecture 22: Deu+eron Quadrupole Moment 24 Nov. 2005 


Assume that the charge distribution is an ellipsoid of revolution with symmetry axis 
along the total angular momentum vector J: 



As for the magnetic dipole moment, 
we specify the intrinsic electric 
quadrupole moment as the expectation 
value when J is maximally aligned with z: 


/ A \ 


Qint = (^2 J 

mj=J 


n m J J 

"Quantum geometry": cos p - yyj - ^ 


If an electric field gradient is applied along the z-axis as shown, the observable 
energy will shift by an amount corresponding to the intrinsic quadrupole moment 
transformed to a coordinate system rotated through angle p to align with the z-axis 



Result: 
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Let Q| ab be the electric quadrupole moment we measure for the ellipsoidal charge 
distribution with mj = J: 



How do we apply this to anything? 

1. A spherically symmetric state (L = 0) has Q jnt = 0 (e.g. deuteron S-state) 

2. Even a distorted state with J = i will not have an observable quadrupole moment 

3. J = 1, L = 2 is the smallest value of total angular momentum for which we can observe 
a nonzero quadrupole moment 

-> these are the quantum numbers for the deuteron D- state! 


Deuteron intrinsic quadrupole moment: 


recall our model for the deuteron wave function: 




a 
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result for the quadrupole moment: 


Q,„, = — \a'b\ ( r 2 ) $D 


note cancellation here 


/ 
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1 b 2 (r 2 ) 

/ DD 
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= + 0.00286 ± 0.00003 bn 


A good model of the N-N interaction can fit both the magnetic moment and the 
quadrupole moment of the deuteron with the same values of a and b! -> 



Basic features of the N-N potential, via the deuteron, etc.: 


1. Independent of the value of L, the state with intrinsic spins coupled to S = 1 
has lower energy 

-> this implies a term proportional to: / has lower energy 

- (^-s,) = - i(s 2 - Sl 2 - s|) = 


2. The deuteron quadrupole moment implies a non-central component, i.e. the potential 
is not spherically symmetric. Since the symmetry axis for Q is along J, Q > 0 
means that the matter distribution is stretched out along the J - axis: 


-1/4, S = 1 
+ 3/4, S=0 




N-N interaction continued.... 


5 


3. There is also a spin-orbit term, as deduced from N-N scattering experiments with 
a polarized beam: 

V s _o - (L.S) 


(This plays a very important role also in determining the correct order of energy 
levels in nuclear spectra - more later!) 


4. Finally, all contributions to the N-N interaction are based on a microscopic meson 
exchange mechanism: 





Where M is a n, p, co ... meson, etc. 

and each term has a spatial dependence 
of the form: 

—m r 

e 

V(r) = g x(spin function ) 

r 



State of the art N-N interaction model: (not to scare anybody...) 6 
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Brace yourself - it takes 2 pages just to write all the terms down!!! 


pseudoscalar mesons: 


V ps( m ps’ r ) = To 


1 £ps 


12 4t r mps 


m 


ps 


m 


Y(m ps r) <r 1 - a 2 + Z(m DS r) S 


ps / 12 




(tensor operator) 


(F.6) 


scalar mesons: 


(they use "a" for spin) 


£ s 


V,(m s ,r) = - ^m s 


l-ip 

4 V m , 


2-| 


1 r-r2i 


y(m/) + 4T 1 [F/y(m » r) + + i z i(»v)*- • s h 


vector mesons: 


/ (F.7) 

(spin-orbit interaction) 


KK> 0 = 7Z 


^ I 

4tt my l 


1+ T ( — 


2 ~\ 


Y(m v r) - — 2 [V 2 Y(m v r) + Y(m v r)V 2 ] 


+ - 


1 l m. 


+ -T 


6 \ m 

1 Sv/v 


) Y(m v r)or • o- 2 - f Z^m^L • S - i Z(m v r)S 12 


2 4tt 


m v {(m v /m) 2 Y(m v r) + |(m v /m) 2 Y(m v r)or • <r 2 


-4Zj(m v r)L • S - 5 Z(m v r),S 12 } 

f 2 

+ 4^ w v {K ^ v /m ) 2 y ( m v r )»'i ‘ ^2 “ AZ(m v r)S 12 } , 



more of this... 
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with 


Y(x) = &~ x /x , 


Z( x) = (mjm) 2 ( 1 + 3/x + 3 lx 2 )Y(x) , 



"Yukawa functions" 
and derivatives... 


Z 1 (x) = 



2 

( 1 lx + 1 lx 1 )Y{x) , 



(F.8) 


and 



We use units such that h = c = 1 (he = 197.3286 MeV fm). The use of the form factor, eq. (3.3), at each 
vertex (with n t( = 1) leads to the following extended expressions: 


K.W = K(m a , r) - 


Ala ~ m l 

A la ~ Ala 


K(4, i,r) + 


A l , i ~ 

Ala ~ <i 


V a (\a > . 


(F.9) 


where A. , = A, + s, = A, — s, e = 10 MeV is an appropriate ch oice 


The full NN potential is the sum of the contributions from six mesons: 
Hr) = 1 V„(r) 

a = / iT,p,T),a),8,a 


\ 



Parameters and Predictions: 
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R. Machleidt et ai, The Bonn meson-exchange model for the nucleon-nucleon interaction 


Table 14 


potential (OBEPR) 

deuteron properties: 

[fjg J 

m a (MeV) A. (GeV) 

LEP 

Theory 

TT 14.9 

138.03 

1.3 

e d (MeV) , 

2.2246 




Pn (%) 

4.81 

p 0.95; [6.1] 

769 

1.3 

2d (fm 2 ) 

0.274 




Md (f%) 

0.8524 

n 3 

548.8 

1.5 

-4 S (fm I;I ) 

0.8860 




D/S 

0.0260 

o) 20; [0.0] 

782.6 

1.5 

G (fm) 

1.9691 




a s (fm) 

-23.751 

5 2.6713 

983 

2.0 

r._ (fm) 

2.662 




a t (fm) 

5.423 

a 7.7823 a 

550 8 

2.0 

r, (fm) 

1.759 

Impressively good agreement for ~ 

10,000 experimental 


low energy scattering 

data points in assorted n- 

-p and p-p scattering 


parameters, etc: 



What does the NN potential look like? 10 


It looks different in different spectroscopic states of the 2N system! 



Only the deuteron is bound! Its quantum numbers have the deepest potential well. 




On to Finite Nuclei! 

Review: 

1. Nuclear isotope chart: (lecture 1) 

• 304 isotopes with ti > 10 9 yrs 

(age of the earth) 

• 177 have even-Z, even-N and J" = 0 + 

• 121 are even-odd and only 6 are odd -odd 

• N » Z for light nuclei and N > Z for 

heavy nuclei 

2. Elastic scattering of electrons: (lecture 7) 

• nuclei are approximately spherical 

1/ 3 

• RMS charge radius R = 1.2 A fm fitted to electron scattering data 

• mass ~ A, and radius ~ A 1/s so the density M/V « constant for nuclei 

(* 2 x 10 17 kg/m 3 ), implying that nuclear matter is like an incompressible fluid 




p ch (r) (fm)~ 


Recall from lecture 7 — Nuclear charge distributions from experiment: 



Fig. 4.3 The electric charge density of three nuclei as fitted by 
p ch (r) = pch/P + exp((r - R)/a)]. The parameters are taken from the compila- 
tion in Barrett, R. C. & Jackson, D. F. (1977), Nuclear Sizes and Structure , 
Oxford: Clarendon Press. 



continued... 
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3, Inelastic electron scattering: (lecture 9) 

• Excited states can be identified, on a scale of a few MeV above the ground state, e.g. 


E (MeV) 


9.641 

1 

3 " 

7.6542. _ 


o + 

4.4389] 

r * 

r 


2 + 








4. Quantum numbers for nuclear states: 


• total angular momentum J, parity n 

• isospin, I: (lecture 13) 

for a nucleus, - \ (Z-N) and I = | m r | , ie lowest energy has smallest I 

• Example: " isobar ic triplet" 14 C, 14 N, 14 0: 


6.09 1 ” 




Nuclear states in general: 15 



• states have integer or half-integer J depending on whether A is even or odd 

• different systematics and energy level spacings for different nuclei 

• some nuclei exhibit "single particle" and others "collective" excitations 

-> different models to describe this complementary behavior 



What is the potential energy function V(r) that nuclei are eigenstates of? 


This is not an easy question! The N-N interaction is too complicated to 
solve in a many-body system: state-of-the-art can go up to A = 3! 

First approximation: a square well potential, width approx, equal to nuclear radius R: 


Assume somehow that we can treat 
the binding of neutrons and protons 
like electrons in atoms - individual 
nucleons have wave functions that are 
eigenstates of some average nuclear 
potential V(r). 

Each nucleon has a binding energy B as 
shown (E = -B) 

Kinetic energy K = (V 0 - B) 

R « 1.2 A 1/3 f m; most of the wave 
function is contained inside the well, 
so this should be approximately the 
right nuclear size... 




Kinetic energy of bound nucleons: 


Key point: once we specify the width of the well, the nucleons are conf ined, and so their 
kinetic energy is essentially determined by the uncertainty principle: 

Confining box of side 2 fm. Ap x Ax~Fi 


> 


Conclusion: the motion is non relativistic; K « 15 MeV 



Simple estimate: 





Now we can specify the potential parameters: *8 


What next? 

• We have a complicated system of A 
nucleons. 


• About half of them are protons, so a 
repulsive (+ve energy) term has to 
be added to the square well to 
account for this (~ few MeV) 

How to connect this model to something 
observable? 


Independent particle model: 

• Assume independent particle motion in some average nuclear potential V(r) as shown. 

• Then we can fill the eigenstates of the potential to maximum occupancy to form a 
nucleus, as is done with electrons in atoms (to 1 st order...) 




Connection to average nuclear properties: 
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• The binding energy of each nucleon, in our model, is a few MeV. 


• The potential energy of a bound nucleon is negative, by ~ 0.3% of its rest mass energy, 
which therefore has to show up as a decrease in its mass. 


• For A nucleons, the total binding energy is: 


A 


A 




z=l 


i=l 


mass of nucleus, M 


The average binding energy per nucleon, B/A, can be determined from mass data 
and used to refine a model for V(r); it ranges systematically from about 1-9 MeV 
as a f unction of mass number for the stable isotopes. 


Reference: F&H ch. 16 



Atomic Mass Units: 
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• By convention, we set the mass of the carbon-12 atom as a standard. 

• Denote atomic masses with a "script" M, measured in atomic mass units, U 

M( 12 C) = 12.0000000000 U (exact!) 1 U = 931.494 MeV (expt) 


Calculation for carbon-12: 


m p = 938.2 MeV 
m n = 939.6 MeV 
m e = 0.511 MeV 


■\ 


► 

j 


6xy mj = 11,269.8 MeV 

i 

12 U = 11,178.0 MeV 

1 


BCeQ = V m, -M = 91.8 MeV 

i 


Binding energy per nucleon in 12 C: B/A = 7.8 MeV; 
Contrast to the deuteron 2 H: B/A = 1.1 MeV 



The famous Binding Energy per Nucleon curve: 
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A 



Implications of the B/A curve: 
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Greater binding energy implies lower mass, greater stability. 

Energy is released when configurations of nucleons change to populate the 

larger B/A region -> nuclear energy generation, e.g. 


Fusion reactions at small A release substantial energy because the B/A curve rises 
faster than a straight line at small A: 


2 H+ 2 H => 3 He + n + 3.27 MeV 
=> 3 H + p + 4.03 MeV 
2 H + 3 H => 4 He + n + 17.6 MeV 






J 


Binding energy 
of products is greater 
than the sum of the 
binding energies 
of the initial species. 


Fission reactions at large A release energy because the products have greater 
binding energy per nucleon than the initial species: 


235 


U + n 


ioo X + 133 Y + 3n + 200 MeV 


distribution of final products 



A semi-empirical model for nuclear binding energies: 
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1. Volume and Surface terms: 

First consider a 1-dimensional row of nucleons with interaction energy per pair = s 


-OO'OOOOO'OO 




V* 

total: A 




A 

i-1 ‘ 

BM— M—M | 

le - A = 2s A - L 

\ 





each has 2 neighbors correction 

for the ends 



B r, A 

— = 2s 

A A 


Approximately constant, with end 
effects relatively smaller at large A. 


By analogy, for a 3-d nucleus, there should be both volume and surface terms 
with the opposite sign, the surface nucleons having less binding energy: 




2. Coulomb term: 
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for a uniform sphere, 


total charge + Ze 

This effect increases the total energy and so decreases the binding energy. 
Simple model: A B = - a c Z 2 A _1/3 




But this is not quite right, because in a sense it includes the Coulomb self energy 
of a single proton by accounting for the integral from 0 to r p ~ 0.8 fm. The nucleus 
has fuzzy edges anyway, so we will have to fit the coeff icient a c to mass data. 


Solution: let AB scale as the number of proton pairs and include a term: 

AB = - a c Z (Z -1) A' 113 => - a c Z (Z -1) A -473 

/l 



3. Symmetry Term: c 

So far, our formula doesn't account for the tendency for light nuclei to have Z = N. 
The nuclear binding energy ultimately results from f illing allowed energy levels in a 
potential well V(r). The most efficient way to fill these levels is with Z = N: 

Simplest model: identical nucleons as a Fermi gas, i.e. noninteracting spin- \ 
particles in a box. Two can occupy each energy level. The level spacing ~ 1/A. 

A mismatch between Z and N costs an energy price of AE at fixed A as shown. 


>" A 

cn 

£_ 

CD 

c 

<D 


O A p 



\r 


t 

At 

A 

ri 

n u 

neutrons 

protons 

states full below 
here.. 




AB = — a A (Z - N) 2 A -1 = -a A (A- 2Z) 2 A' 1 



4. Pairing Term: 


Finally, recall from slide 1 that for the case of even A, there are 177 stable 
nuclei with Z and N both even, and only 6 with Z and N both odd. Why? 

-> Conf igurations for which protons and neutrons separately can form pairs 
must be much more stable. All the even-even cases have J* = 0 + , implying 
that neutrons and protons have lower energy when paired to total angular 
momentum zero. 

Solution: add an empirical pairing term to the binding energy formula: 


A B nnir = 8 = 

~+l~ 

0 

a A~ 3/4 

pair 

_-l_ 

P 


with +1 for even-even, 0 for even-odd, and -1 for odd-odd 

Full expression: 

B(Z,A ) = a v A-a s A 2/3 -a c Z (Z-l) A~ 113 - a A (A-2Z) 2 A -1 + 8 



Fitting of coeff icients to data: 27 


B(Z,A) = a v A-a s A 2/3 -a c Z(Z-1)A~ 1/3 - a A (A—2Z) 2 A -1 + 8 




One more look at the Binding Energy per Nucleon curve: 
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